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We show that the "perfect fluid" elliptic flow of the bulk hadrons and the unexpectedly strong 
quenching of heavy quark jet fragments in Au+Au reactions at 200 AGeV can be simultaneously 
accounted for within leading order AdS/CFT holography with a common large t'Hooft coupling 
A = 9ymN c ~ 30. In contrast, weakly coupled quark-gluon plasma models have so far failed 
to describe the observed correlation between these soft and hard observables even for couplings 
extrapolated to a a ~ 0.5. We show that phenomenological applications of the classical trailing 
string AdS/CFT holographic solution are furthermore remarkably robust to higher order curvature 
corrections C(l/A 3 '' 2 ) in type IIB supergravity theories, as well as to worldsheet fluctuation cor- 
rections (^(l/A 1 / 2 ) that were not considered previously. We emphasize the importance of future 
measurements at RHIC and LHC of the correlation between identified charm and beauty quark 
hard (pr > 10 — 30 GeV) jet quenching observables and low transverse momenta (pr < 1 GeV) 
bulk elliptic flow observables to further tests the limits of applicability of conformal holography to 
strongly coupled quark-gluon plasma. 

PACS numbers: 25.75.-q, 11.25.Tq, 13.87.-a 
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I. INTRODUCTION 

Two of the most remarkable experimental discoveries 
at the Relativistic Heavy Ion Collider (RHIC) are the 
observations of a factor five quenching of (high trans- 
verse momentum or high quark mass) jets and the nearly 
"perfect fluid" elliptic flow of low transverse momentum 
hadrons [l|. These observations and bulk multiplicity 
production systematics have been interpreted as provid- 
ing evidence for the formation of two new forms of QCD 
matter 1) a locally equilibrated strongly-coupled Quark- 
Gluon Plasma (sQGP) and 2) its non-Abelian classi- 
cal Color Glass Condensate (CGC) initial nuclear field 
source @, Q. 

However, it has been a challenge to find a single con- 
sistent theoretical framework to explain simultaneously 
both soft (long wavelength) and hard (short wavelength) 
properties of the sQGP phenomena. Attempts to explain 
bulk collective flow 041] based on perturbative QCD 
parton transport approaches 0413 1 require large cou- 
pling extrapolations a s = g\ M j / ^K — > 0.6, which on the 
other hand overestimate of the opacity of the sQGP to 
high transverse momentum jets. At moderate coupling 
a s ~ 0.4 the observed opacity for light quark jets can be 
well accounted for but at the expense of too high viscosity 
and underestimating the elliptic flow. The heavy quark 
jet quenching data pose an especially difficult challenge 
for perturbative QCD energy loss models 

For large values of the QCD coupling, the t'Hooft pa- 
rameter A = gqcD-Nc > 20 (with N c = 3) may already be 
large enough to validate string theory inspired classical 
AdS/CFT holographic models [H [D| of the sQGP. In 
this approximation, the viscosity to entropy ratio is nat- 
urally small [HI and close to its lower unitarity bound 
@; EH ■ Classical strings in the AdS background further- 



more provide a dual holographic model of heavy quark 
drag with dE/dx oc \f\ET 2 /Mq that can easily account 
for the high opacity of the sQGP to heavy quark jets 
[l5j |. The key question that we address in this paper is 
whether there exists a single value of A that could ac- 
count for both the bulk elliptic flow as well as the strong 
heavy quark quenching simultaneously. We test at the 
same time the consistency of the CGC initial field source 
geometry out of which the sQGP forms. Our conclusion 
is that conformal holography can indeed account for the 
observed hard/soft correlation data with the CGC ini- 
tial geometry with A ~ 30 while classical Glauber initial 
geometry is not consistent within present experimental 
errors. We emphasize the importance of testing future 
identified charm and bottom jet quenching data system- 
atics and their correlation with soft bulk flow observables. 



II. STRONGLY-COUPLED Af = 4 
SUPERSYMMETRIC YANG-MILLS AND 
HOLOGRAPHY 

The planar limit of strongly-coupled M = 4 Supersym- 
metric Yang-Mills (SYM) theory provides a powerful phe- 
nomenological model of the sQGP at zero baryonic chem- 
ical potential in the temperature range where the trace 
anomaly is small. Soft collective flow phenomena are con- 
trolled by the system's entropy density, s(T), and shear 
viscosity, r](T). In this conformal c 2 . = dhiT/dhis = 1/3 
theory, the bulk viscosity vanishes and the dimensionless 
ratio rj/s is independent of temperature. While confor- 
mal invariance is broken in real world QCD, the confor- 
mal symmetry of M = 4 SYM played an important role 
in the conjectured classical supergravity dual description 
[l2j of this CFT on a AdSs x S5 curved spacetime. In 
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the N c , A — > oo limits, the compactification of a lOd type 
IIB string theory on a 5d sphere of radius L leads to a 
5d Einstein action with a negative cosmological constant 
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The effective 5d gravitational coupling is taken as G5 ~ 
1/-/V 2 <C 1, and the AdS's metric G M „ is found to be a 
stationary solution with curvature 1Z = —12/L 2 . The 
t'Hooft coupling in the gauge theory is identified with 
L 2 /a', where Va* — £ s is the fundamental lOd string 
length. The a' expansion in the gravity dual description 
is mapped into a series in 1/y/X in the gauge theory [l3j . 
The extra • • • terms in Eq. (jXJ) correspond to the correc- 
tions to the leading order supergravity action that appear 
when L 2 /a' and/or N c are large but finite [13j . The lead- 
ing "stringy" corrections to type IIB string theory were 
found to come from terms ~ o/ 3 7?. 4 to the tree level ef- 
fective action [13j and they induce contributions of order 
1 / A 3 / 2 to the quantities mentioned above [TH [l|| [l?} ■ 

Breaking of supersymmetry on the gauge theory side 
can lead in the gravity dual description to lower or- 
der quadratic or cubic curvature corrections. A class 
of Gauss-Bonnet generalizations of the effective 5d Ein- 
stein action was considered in fl8l-20) . These corrections 
are characterized by a second assumed small dimension- 
less parameter, Xqb ~ 1/N C , which is related to the 
central charges c and a that characterize the conformal 
anomaly in curved spacetime of the dual CFT as noted 
in Eq. (2.14) of Ref. [l8j]. Varying Xqb provides a para- 
metric way to explore small deformations of the original 
Af = 4 SU(N C ) SYM theory. Interest in possible Gauss- 
Bonnet deformations of conformal holography arose when 
Kats and Petrov [H found that for Af = 2 Sp(AU 
SYM, X GB = 1/8N C , and the KSS viscosity bound (14) . 
r l/ s > 1/4tt, could be violated by 17% for N c = 3, in 
such models. Ref. UH further showed that a large class 
of other effective CFTs lead to similar Agb oc l/N c ef- 
fects. In Refs. [HIHH it was found however that causality 
and positive energy flow limit deformations to a narrow 
parametric range —7/36 < Xqb < 9/100. We find below 
that RHIC hard/soft correlation observables are in fact 
compatible with Xqb = 0. 

Conformal holography refers to the predicted temper- 
ature independence of s/ssb, the entropy density of a 
very strongly coupled SYM compared to its ideal Ste- 
fan Boltzmann limit, and also of the viscosity to entropy 
density ratio. The heavy quark jet relaxation rate, 1/tq, 
is controlled by fiQ = \fX~TiT 2 /2Mq for a heavy quark 
with mass Mq in a plasma of temperature T [15|, [22| . 
The relaxation time is related to the heavy quark energy 
loss per unit length through tq(X) = —l/(dlogp/dt) = 
—l/(dlogE/dx), where p = Mqjv and v = p/E. 

Our analysis is based on the following remarkably 
simple algebraic expressions relating three fundamental 
properties of large N c , Af — 4 conformal SYM plasmas 



at large t'Hooft coupling A: 
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where C3 = 15C(3) sa 18, K\ ~ — 1 is a new worldsheet 
fluctuation correction amplitude discussed below. The 
A~ 3 / 2 correction to the entropy density ratio in Eq. @ 
was found in [l3T | while the analogous correction to 77/s 
was obtained from [l^. The finite t'Hooft coupling cor- 
rection to the heavy quark energy loss, C3A 3 / 2 /16, is a 
new result reported here (see the derivation below) and 
it is needed for a consistent application of the strongly- 
coupled M = 4 SYM model to heavy ion reactions, as 
will be argued below. 



III. HEAVY QUARK JETS IN Af = 4 SYM AT 
LARGE (BUT FINITE) T'HOOFT COUPLING 

The hard observable we consider here is heavy quark 
energy loss. An object of mass Mq in the fundamen- 
tal representation of SU(N C ) added to SYM is described 
holographically using a D7 brane [23J that ends at a ra- 
dial direction u m ~ Mq away from the black brane hori- 
zon Uh ~ T '. A heavy quark moving with constant ve- 
locity v is then described in terms of a curved string in 
the radial direction which has the upper end attached to 
the bottom of a D7 brane while the lower end remains 
attached to the black brane [lg, [22J . When N c — > 00 and 
A > 1 (but finite) the string dynamics is embedded in a 
black brane curved background spacetime of the type 



ds 2 — Goo(u)dt 2 + G xx (u)dx 2 + G uu (u)du 2 (5) 



where 
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G uu {u) = ^ (l + 0(c/ 3 )) 
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and G xx — u 2 / L 2 (l + 0(a' 3 )) (the finite a' corrections 
to the metric we use can be explicitly found in |13l|). 

The classical Nambu-Goto string action (with dilaton 
ommitted) is 
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where 



g = dot g ab = G^daX^dbX" 
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is the induced worldsheet metric, a a = (r, a) are the in- 
ternal worldsheet coordinates, G lll/ {X) is the background 
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metric, and X^ — X^(t, a) is the embedding of the string 
in spacetime. By restricting the analysis to the classical 
string the energy loss always becomes proportional to yA 
at leading order because of the string tension l/(27ra') 
in the NG action. 

The trailing string ansatz (where r = t, a — u and 
X^it, u) = (t, x + vt + £(u), 0, 0, u)) is assumed to de- 
scribe the asymptotic behavior of a string whose endpoint 
("the heavy quark") moves with velocity v in the x di- 
rection and located at a fixed AdS radial coordinate near 
the boundary u m ^> uh [HI, HH . The black brane hori- 
zon coordinate uh oc Ta' is determined by Goo(uh) = 0. 
Using the ansatz above and the string's classical equa- 
tions of motion, it can be shown [15[ that the drag force 
dp/dt = —Cv/(2ira'), where C is a constant determined 
by the negativity condition that 



g(u) = G uu (Goo + v 2 G xx ) (l + ^ ) < (10) 



for Uh < u < u m . However, both the numerator and de- 
nominator in Eg. (|10l) change their sign simultaneously 
at a certain u* [15| given by the root of the equation 
G o(u*) + v 2 G xx (u*) = 0. This fixes C = G xx (u*) 
and dp/dt = — v G xx {u*) / {2-koi'). Neglecting higher- 
order derivative corrections in J\f = 4 SYM one finds 
u* = Uhy/j, where 7 = l/\/l — v 2 . The condition that 
u* < u m leads to a maximum "speed limit" for the heavy 
quark jet to be consistent with this trailing string ansatz 
given by 7max < u 2 n /u 2 h 0]. 

Using the metric derived in [l3] to 0(a' 3 ) and 
the classical NG action, one can compute the effects 
of quartic corrections on the classical drag force and 
determine it* perturbatively to 0(A~ 3 / 2 ) as u* = 
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and the drag force 
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The analytical result for the classical heavy quark drag 
force at large (but finite) t'Hooft coupling in J\f = 4 SYM 
computed above is one of the main results of this paper. 
The heavy quark mass at T = is Mq = u m j \2ira') 

and, to leading order in 1/A, u 2 m ju\ ~ -^r- Thus, the 
corrected u* displayed above defines a new speed limit 

7m ~ [l — jg (^7^ — l)] 1 after neglecting terms of 
0(1/7, l/iV c ). Note that j m and dp/dt decrease with 
increasing 77/ s. 

The drag force experienced by a heavy quark in the 
strongly-coupled J\f = 4 SYM plasma can also receive fi- 
nite A corrections from quantum fluctuations of the string 
worldsheet, as it will be discussed below. 



A. Finite t'Hooft coupling corrections to heavy 
quark energy loss 

While the large amount of supersymmetries in M = 4 
SYM excludes the corrections to r]/s due to 1Z 2 and 7Z 3 - 
like terms in its lOd gravity dual action 13(, this is not 
the case for observables such as the string drag or the 
effective heavy quark potential (25l - [27| . There are no 
known reasons to exclude terms in the a' expansion of the 
2d world sheet theory. Physically, these corrections corre- 
spond to fluctuations of the string worldsheet around its 
minimum area. For instance, the a' corrections just from 
the worldsheet string loops defined on top of a supergrav- 
ity background (where the gravity dual's action contain 
only quadratic terms in the spacetime derivatives) would 
contribute to the energy loss E of a heavy quark in the 
medium as follows 
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where x = (irT)x is a dimensionless quantity, E = Mq7, 
and Ki with i = 1 , . . . , 4 denote the contribution from 
i — th string loop corrections. The first term outside 
the parenthesis is the standard result obtained before in 
the A — > 00 limit [IH, [22[. For physical quark masses 
the classical string description is only applicable when 

-t? T7— < 1, which indicates that A cannot indeed be 

2 Mq 

taken to be infinitely large. In a critical string theory 
such as type IIB, the loop coefficients are expected to be 
finite and, in general, v dependent. The 1-loop coefficient 
of the correction to the heavy quark potential in the vac- 
uum (which gives a term analogous to the Ki/vA above) 
was recently found numerically to be ~ —1.33 [27l ]. 

As shown in Eq. (fTTj) . the a' 3 correction to the back- 
ground studied in [13| affects the energy loss. This con- 
tribution can be taken into account independently of the 
worldsheet fluctuations at least to lowest order. The ul- 
trarelativistic limit of Eq. (fTTj) is 
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This is the new term that was included in Eq. (j4j . Note 
that with this type of correction suggests that the mag- 
nitude of the energy loss becomes larger than the value 
obtained at infinite coupling as A decreases (or rj/s in- 
creases), which is certainly counterintuitive and most 
likely unphysical. This poses no harm if after the string 
loop corrections in Eq. (112[) are included the correct phys- 
ical behavior (i.e., smaller coupling means smaller energy 
loss) is recovered. This implies that worldsheet fluctua- 
tions must be taken into account in the study of finite A 
corrections to the energy loss in JV = 4 SYM. 

If one includes both the string loop corrections and the 
bulk-induced new term one obtains (assuming that the 
v — > 1 limit of the loop corrections is well defined) 
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A comparison between Eqs. fl2J), fl3J), and ([T4"]l shows that 
the heavy quark energy loss can be in principle much 
more sensitive to finite A corrections than the bulk quan- 
tities. This is also true, for instance, for heavy quark 
bound states described within AdS/CFT On the 

other hand, bulk quantities such as rj/s are finite at any 
(large) value of the coupling. 

For our applications, we consider the range A ~ 5 — 30 
and N c = 3. Also, we neglect other formally higher order 
corrections [l7[ ~ 0(V\/N 2 ) to the Af = 4 SYM entropy 
density and shear viscosity and set n 2 = K3 = K4 = 
in this first attempt to test conformal holography of 
hard and soft observable correlations in high energy A+A 
collisions. 



IV. CORRELATION BETWEEN THE HEAVY 
QUARK JET Raa AND BULK ELLIPTIC FLOW 

The dependence of soft phenomena on rj/ s in nearly 
perfect fluids can be studied using viscous relativistic hy- 
drodynamics [29j . Once v 2 as a function of r//s is known 
from hydrodynamic simulations, one can use Eq. ([3]) to 
obtain its dependence on A and estimate how much ellip- 
tic flow is generated in a strongly-coupled SYM plasma. 
In order to compute v 2 (pt, v/ s ) f° r a given centrality in- 
terval C, we employ a linear fit to the numerical results 
of for both Glauber [3(3] and CGC initial transverse 

profiles mm 

wa(PT)f?/s,C) = a(p T ) e 2 (C) (1 - a x rj/s) (15) 



where e 2 (C) = (y 2 



y 2 )c 



is the average 



initial elliptic geometric eccentricity for the centrality 
class C. We consider the 20-60% centrality class be- 
cause, as shown in Fig. 23 of there is good agree- 
ment at p T ~ 1 GeV between STAR u 2 (4) and PHENIX 
v 2 (BBC) data and non-flow effects [33| are reduced. 
To rescale the minimum bias results of Ref. [29| to 
the considered 20-60% centrality class we use the factor 



-Glaub 



(20-60%)/ 



-Glaub 



(0 - 92%) = 0.317/0.281 = 1.128 



from Ref. [6j. Our fit to the rescaled numerical results 
of HI gives ax ~ 2.5 and a(p T = l)e 2 (20 - 60%) « 
0.14 (0.098) for CGC (Glauber) initial conditions. In 
this case, the Glauber initial conditions do not generate 
enough elliptic flow and only the CGC initial conditions 
can fully describe the elliptic flow data. 

The nuclear modification factor of single non-photonic 
electrons, R AA (pr), obtained from quenched heavy quark 
jets can be computed at strong coupling [34| by using the 
generalized energy loss in Eq. (fi"4|) to compute the path 
length dependent heavy quark fractional energy loss e. In 



this case Raa = ((1 - e)" Q 
dependent spectral index uq 



, where uq(pt) is the flavor 
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obtained from FONLL production cross sections [35| . 
The path length average of the nuclear modification at 
impact parameter b is computed using a Woods-Saxon 
nuclear density profile with Glauber profiles Ta(x±) with 
&NN = 42 mb. For 0-10% centrality triggered data both 
Glauber and CGC geometries lead to similar numerical 
results for Raa [3l| . The distribution of initial hard jet 
production points at a given x± and azimuthal direction 
<j> is taken to be proportional to the binary parton colli- 
sion density, Taa(x±, b). 

We consider a longitudinally expanding local (partic- 
ipant) parton density p(x±,b) = xPpart(x±,b)/r, where 
X = (dNir I dy) I Npmt and p par t is the Glauber partici- 
pant nucleon profile density. In our calculations we use 
a reduced temperature T C ft = 0.74(5/ Ssb) 1/3 T QC d to 
take into account the fewer number of degrees of free- 
dom in a strongly-coupled QCD plasma, which is similar 
to the prescription given in [36j. The heavy quark mod- 
ification factor is 

/•27T 



RY*(pr,b) = I d(f> I d 2 x 



T A A(x±,b) 
■ 2^Y bin (6) 



x exp 
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where Yb; n is the number of binary collisions and the 
relaxation time tq is defined in terms of (|14l) . 

Here, To = 1 fm/c is the assumed plasma equilibra- 
tion time and Tf is determined from T(£, r/) = Tf = 140 
MeV, i.e, the time at which the local temperature falls 
below a freeze-out temperature taken from [29(. We as- 
sume that the non-photonic electron modification factor 



is R A a(pt 
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where 



Pt/0.7 is an estimate of fragmentation effects fllj j . 
The Af = 4 SYM prediction for the soft v 2 and hard 
Raa correlation, computed to lowest order in the cou- 
pling (supergravity limit), is shown in Fig. 1 (horizontal, 
red line) for A = 5, . . . , 30. Only the results obtained 
with CGC initial conditions are shown. The eccentricity 
of the Glauber initial geometry models are too low to fit 
the data. In this case inrj/s — 1 and the energy loss is 
~ VX [HI • The agreement between this lowest order 
calculation and the RHIC data is rather remarkable for 
large A ~ 20 — 30. While finite coupling corrections to 
r]/s only slightly shift v 2 to the left, R%a involves an ex- 
ponential of the energy loss and, thus, one may expect 
this quantity to display a stronger sensitivity to 1/A" cor- 
rections. While most of the coefficients in the expansion 
(fl4l are unknown, a useful check of the robustness of the 
lowest order result can be done by repeating the calcu- 
lations with the known 1/A 3 / 2 corrections to rj/s in Eq. 
© and approximating Eq. (|14p with the 1-loop world- 
sheet correction coefficient to the heavy quark potential 
in vacuum Ki = —4/3 found in [27] (with other correc- 
tions neglected). This leads to the dashed blue curve 
in Fig. 1. Note that because A ~ 30 is so large even 
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FIG. 1. (Color online) The observed correlation between the 
nuclear modification, Raa(Pt = 5.5), of 0-10% centrality non- 
photonic electrons from heavy quark jets and the bulk elliptic 
flow moment V2{pt = 1 GeV) of pions produced in 20-60% 
semi-central reactions is shown by the orange (grey) ellipse 
from PHENIX (STAR) measurements of Au+Au 200 AGeV 
data at RHIC [j, 133 . 140] ]. Curves show predictions based 
on conformal holography with t'Hooft coupling in the range 
5 < A < 30 and assuming the initial participant geometry 
predicted by the KLN Color Glass Model HlMlH and with 
total entropy constrained by the measured total pion multi- 
plicity dN^/dy = 1000 in central collisions. The horizontal 
black line corresponds to the lowest order supergravity result 
with [k\ = 0,C3 = 0) in Eqs. (2-4). The green curve is for 
(kj = 0, C3, ~ 18) including only 1Z 4 stringy curvature cor- 
rections to the TZ 1 Einstein action. The dashed blue curve 
includes, in addition, an estimate for the leading worldsheet 
fluctuation correction W 1 with K\ = —1.33 (see text). The 
purple band denotes a range of these observables far from the 
observed data as computed in pQCD-based transport models 



this formally much larger correction is negligible com- 
pared to the leading supergravity prediction. We have 
checked [4l[ that adding phenomenologically a "pertur- 
bative" quadratic Gauss-Bonnet curvature corrections to 
the Einstein action dual to TV = 4 SYM (T^|20j including 
also its modification of the heavy quark energy loss from 
[13, HI], following the suggestion made in Ref. [18], does 
not lead to any significant changes in the results shown 
in Fig. 1 for the range |Agb| < 0.09 allowed by causality 
and positivity of energy flux. This indicates that possi- 
ble 1/N C corrections do not alter the conclusions of this 
paper. 

The light purple band near the top right corner of 
Fig. 1 shows the range of hard/soft correlations that far 
from the data predicted by pQCD based transport mod- 
els [1, |Tl|. Our results suggest that the sQGP discovered 
at RHIC are remarkably well described at both large and 



small wavelengths by an approximate conformal holo- 
graphic model, while the observed hard/soft correlation 
severely challenges perturbative QCD approximations. 

We conclude by demonstrating that future identified 
charm and beauty jet quenching data at RHIC and LHC 
will soon allow much more stringent and sensitive test of 
holographic QCD phenomenology. The predicted small 
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FIG. 2. (Color online) The predicted correlation between 
identified heavy quark jet quenching, R^ A , and soft elliptic 
flow, V2(pr = 1 GeV), based on approximate R 1 + W 1 + R 4 
conformal holography that includes leading order worldsheet 
fluctuations and quartic curvature corrections via Eqs. (2- 
4). The upper red curve corresponds to bottom quark flavor 
tagged jets while the lower green curve corresponds to charm 
flavor tagged jets. The blue dashed curve is the nuclear mod- 
ification of single electrons from b and c jets estimated from 
R*aa ~ 0.6R b A A + 0AR c aa with p? m 0.7p^, as also shown in 
Fig. 1. 

double ratio R AA /R b AA <§; 1 is a robus t sig nature of AdS 
string drag models as emphasized in [34]], which differs 
significantly from pQCD-based transport models and can 
be readily tested once flavor tagged jet measurement in 
A+A reactions become feasible at RHIC and LHC. 

An important open theoretical problem is to generalize 
the above analysis to non-conformal holographic models 
that that take into account the conformal anomaly near 
T c ~ 170 MeV as predicted by lattice QCD. Furthermore 
a more consistent and quantitative light quark/gluon 
jet non-conformal holographic QCD phenomenology will 
need to be developed to account simulaneously also for 
the observed nuclear modification high px pions from 
light quark and gluon jets, R aa (j>t)- The holographic 
string drag model postulated for heavy quarks is inap- 
plicable for light quark jets and may break down already 
for charm quark jets. This possibility underlines the need 
for the next generation experiments at RHIC and LHC 
that can measure simultaneously charm and bottom jet 
quenching observables as well as pion nuclear modifica- 
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tion factors and their correlations with the bulk elliptic 
flow and total entropy production. 
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